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ABSTRACT 

We present a study of the fundamental plane, FP, for a sample of 7 1 dwarf 
galaxies in the core of Coma cluster in magnitude range —21 < Mj < — 15. 
Taking advantage of high resolution DEIMOS spectrograph on Keck II for mea- 
suring the internal velocity dispersion of galaxies and high resolution imaging of 
HST/ACS, which allows an accurate surface brightness modeling, we extend the 
fundamental plane (FP) of galaxies to ~1 magnitude fainter luminosities than all 
the previous studies of the FP in Coma cluster. We find that, the scatter about the 
FP depends on the faint-end luminosity cutoff, such that the scatter increases for 
fainter galaxies. The residual from the FP correlates with the galaxy colour, with 
bluer galaxies showing larger residuals from FP. 

We find M/L oc M-° 15±0 - 22 in F814W-band indicating that in faint dwarf 
ellipticals, the M/L ratio is insensitive to the mass. We find that less massive 
dwarf ellipticals are bluer than their brighter counterparts, possibly indicating 
ongoing star formation activity. Although tidal encounters and harassment can 
play a part in removing stars and dark matter from the galaxy, we believe that the 
dominant effect will be the stellar wind associated with the star formation, which 
will remove material from the galaxy resulting in larger M/L ratios. We attribute 
the deviation of a number of faint blue dwarfs from the FP of brighter ellipticals 
to this effect. 

We also study other scaling relations involving galaxy photometric properties 
including the photometric plane. We show that, compared to the FP, the scatter 
about the photometric plane is smaller at the faint end. 

Key words: galaxies: clusters: individual: Coma; galaxies: elliptical and lenticu- 
lar, cD; galaxies: dwarf; galaxies: kinematics and dynamics; galaxies: fundamen- 
tal parameters; galaxies: evolution 



* Based in part on observations made with the NASA/ESA Hub- 
ble Space Telescope, obtained from the data archive at the Space 
Telescope Institute. STScI is operated by the association of Uni- 
versities for Research in Astronomy, Inc. under the NASA contract 
NAS 5-26555. These observations are associated with programme 
GO10861. 

f Some of the data presented herein were obtained at the W.M. 
Keck Observatory, which is operated as a scientific partnership 
among the California Institute of Technology, the University of Cal- 
ifornia and the National Aeronautics and Space Administration. The 



1 INTRODUCTION 

The photometric and dynamical properties of the elliptical 
galaxies could be explained by multivariate analysis which 
reveals the fundamental scaling relations of these galaxies. 
One of the most well-known manifolds showing a tight cor- 
relation between R e (i.e. the radius encompassing half-light 
of the galaxy), I e (i.e. the mean surface brightness within R e 
in flux unit), and a (i.e. the galaxy internal velocity disper- 

Observatory was made possible by the generous financial support of 
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Figure 1. Comparing the Sersic parameters in two different filters 
(F814W & F475W). For both filters, the initial parameters to run 
Galfit are chosen from the SExtractor catalogue of F8f4W images. 
The solid black line is the locus of ngi4 = 71475. The blue dashed 
line is the best linear regression between the two parameters. 71814 
values are 10% higher than 7*475 values, on average. Red triangles 
and green diamonds represent galaxies of Co09 and MG05 samples, 
respectively. Black asterisks are dwarf galaxies of our DEIMOS 
sample. 



sion), is the Fundamental Plane (hereafter FP; Djorgovski & 
Davis 1987; Faber et al. 1987, Dressier et al. 1987; Bernardi 
et al. 2003). 

In theory, the FP is derived from the virial theorem as 
R e oc a 1 'Ie 1 (M 1 'L)" 1 , where I e is the effective surface 
brightness in flux units, calculated within the half-light ra- 
dius, R e , of the galaxy, a is the galaxy internal velocity dis- 
persion and M/L is its mass-to-light ratio. Assuming that 
the M/L is expressed by a power-law function of a and/or 
the effective surface brightness, I e , the FP relation is simpli- 
fied as 

log(R e ) = A log{a) + B (/x) e + C, (1) 

where (/i) e is the mean surface brightness in 
mag / arcsec 2 unit and is defined as (/j) e = —2.5log(I e ) + 
cte. Although the shape of the FP and its coefficients dif- 
fers for different gravitationally bound systems from globu- 
lar clusters (Burstein et al. 1997) to galaxy clusters (Schaef- 
fer et al. 1993; Fritsch & Buchert 1999; Zaritsky et al. 2006a: 
ZGZ06), there is no doubt about its existence (Lucey, Bower 
& Ellis 1991 and its references). In reality, the coefficients 
of the FP relation differ from the prediction of the virial 
theorem. The observed coefficients are A=1.24 and B=0.33 
(J0rgensen, Franks and Kjasrgaard 1996: JFK96) while the 
virial theorem predicts A=2.0 and B=0.4. This difference, 
often referred to as the "tilt" of the FP, is mainly attributed 
to different formation histories and evolutionary processes. 

The difference in FP coefficients of different spheroidal 
systems with different mass, size and luminosities, can be 
explained by evolution of the M/L ratio as a function 



of stellar population [age, metalicity or initial mass func- 
tion (IMF)] and/or dark matter content (Tortora et al. 2009; 
Grillo & Gobat 2010; Graves & Faber 2010). In addition, 
the absence of homology, i.e, the fact that, the structure 
of spheroids is scalable regardless of their size, can be the 
source of the FP tilt (D'Onofrio et al. 2008; Trujillo, Burk- 
ert & Bell 2004). On the other hand, some authors stud- 
ied the role of dissipation in explaining the nature of the 
FP (Ribeiro & Dantas 2010; Hopkins, Thorns & Hernquist 
2008: HCH08). HCH08 claimed that the non-homology or 
change in the dark matter distribution are not the main 
drivers of FP tilt. Studying the early-type galaxies in 59 
nearby galaxy clusters, D'Onofrio et al. (2008) have found a 
strong correlation between the FP coefficients and the local 
cluster environment and no strong correlations with inter- 
nal galaxy properties (e.g. S'ersic index and galaxy colour). 
Moreover, FP coefficients are independent of global prop- 
erties of clusters such as radius, X-ray luminosity and clus- 
ter velocity dispersion (D'Onofrio et al., 2008). On the other 
hand, while Reda, Forbes & Hau (2005) have shown that iso- 
lated early-type galaxies lie on the same fundamental plane 
as galaxies in high-density environments, cluster galaxies 
have also less intrinsic scatter in their properties compared 
to field galaxies (de Carvalho & Djorgovski, 1992). 

The study of dark matter in dwarf galaxies showed that 
below the critical virial velocity, which is estimated to be 
~100 km s , interstellar gas removal via supernova explo- 
sions become important (Dekel & Silk 1986: DS86). This 
mechanism has been invoked also to explain the shape of the 
low mass dwarf galaxies and their mass profiles (Sanchez- 
Janssen et al. 2010; Governato et al. 2010). Differences be- 
tween the FP of giant and dwarf galaxies have been known 
for some time (Nieto et al. 1990; Bender et al. 1992; Guz- 
man et al. 1993). Peterson & Caldwell (1993), studying a 
sample of nucleated dwarfs, they found a change in M/L 
ratio with luminosity as predicted by the scaling relations of 
DS86. The study of 17 Virgo dwarfs (-17.5< M v <-15.48) 
by Geha et al. (2002 & 2003) placed dEs on a plane parallel 
to, but offset, from that occupied by normal elliptical galax- 
ies. On the other hand, Graham & Guzman (2003), based 
on a detailed surface photometry from HST archival images 
for a sample of dEs in the Coma cluster, claimed that the 
E-dE dichotomy is a consequence of a continuously vary- 
ing profile shape with galaxy luminosity. They attributed the 
E-dE dichotomy to the fact that the light profile of dwarf 
galaxies do not follow the de Vaucouleurs law (1948) and 
therefore, concluded that giant and dwarf ellipticals (here- 
after dEs) were not formed by different mechanisms. How- 
ever, as found by Geha et al. (2003), modelling the dEs with 
Sersic functions (1968) does not explain why dEs lie on a 
different region of the FP. Although it is clear that dwarf el- 
lipticals do not lie on the fundamental plane as defined by 
brighter ellipticals (Geha et al. 2003), Zaritsky et al. (2006 
a,b) showed that the differences between the "fundamen- 
tal planes", defined by different families of spheroids from 
dwarf ellipticals to galaxy clusters, is the result of a non- 
linear relationship between logcr and log M/L. Desroches 
et al. (2007) also found the sensitivity of the coefficient A 
of the FP to the luminosity faint-end cutoff. D'Onofrio et 
al. (2008) were among the first to question the universality 
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Figure 2. Two examples of the fitted Sersic function to the galaxy light profile. Red open triangles show the real profile of the galaxies 
while the blue dotted and black solid curves represent the best fit from a single Sersic and Sersic+gaussian models, respectively. The vertical 
black double-dotted lines show the location of the effective radius, R e , of these galaxies. The horizontal axes show the average radius of the 
ellipsoidal isophotes (r = %/ab). The red dash-dotted curves represent the typical profile of the PSF. Left) It shows the surface brightness 
profile of GMP3080. The central excess light forced us to add an extra gaussian component which results in the best fit with the Sersic index of 
n=1.0 (the black solid curve). The FWHM of the added central Gaussian function is 0.13". The Sersic index corresponding to the blue dotted 
curve is 0.9. Right) The surface brightness profile of GMP3133, which is accurately described with a single Sersic profile. 



of FP and proposed the bending of FP of early-type galax- 
ies. On the other hand, Fraix-Burnet et al. (2010) studied a 
sample of galaxies with the redshift range 0.007-0.053, and 
found that the global FP is a mixture of several fundamental 
planes, with different thickness and orientations, and con- 
cluded that the FP is not a bent surface. They found seven 
groups of fundamental planes which are assigned to differ- 
ent galaxy assembly histories and formation scenarios. 

Another fundamental scaling relation, based on pure 
photometric observables, was introduced by Khosroshahi et 
al (2000 a,b), called the Photometric Plane (hereafter PHP) 
of the galaxies. The Sersic index in the PHP, which replaces 
the velocity dispersion in FP, carries information on the light 
profile in galaxies and it is correlated with the dynamical 
properties of the galaxies such as the total mass and internal 
stellar velocity dispersion (Graham 2002: GR02). Steeper 
the light profile at the centre of the galaxy (higher Sersic in- 
dex), more massive the galaxy and larger the central velocity 
dispersion. 

FP and PHP could be potentially used to derive distance 
to any single galaxy. GR02 has compared both photomet- 
ric and fundamental planes for exactly the same sample of 
early-type galaxies in Virgo and Fornax clusters, and con- 
cluded that the scatter about the PHP is ~ 15-30% more than 
the scatter about FP. Since the PHP is purely based on pho- 
tometric quantities, it can be replaced with FP wherever the 
expensive kinematic data are not available and/or high level 
of accuracy is not required. 

In Paper I in this series (Kourkchi et al. 2011a) we 
used the high spectral resolution of DEIMOS (Faber et al. 
2003) on the Keck II telescope to measure precise veloc- 
ity dispersions for a sample of 41 faint elliptical galaxies in 
Coma cluster. Of this sample, 28 galaxies lie in the observed 
footprints of the HST/ACS Coma Treasury survey (Carter 
et al. 2008). To this sample we add a further 43 galaxies 
with velocity dispersions from Matkovic & Guzman (2005: 
MG05) or Cody et al. (2009: Co09), which also lie in these 
footprints, giving us a combined sample of 71 objects with 
— 22 < Mr < —15 with precise velocity dispersions mea- 



sured from one of these sources. We noticed that fainter 
dwarfs show a departure from brighter ellipticals on the 
Faber Jackson relation (Faber & Jackson 1976) which in- 
dicates they have higher velocity dispersion and therefore 
are more massive or less luminous than is predicted by L-cr 
linear trend. 

Any study of the differences in the FP of dwarf and 
giant ellipticals helps us to test the galaxy formation sce- 
narios. In this paper, we present the structural parameters 
(effective radius, effective surface brightness, central sur- 
face brightness, Sersic index and concentration parameter) 
for galaxies in our sample, based on the Treasury HST/ACS 
images in the F814W and F475W filters. The precise mea- 
sured structural parameters from the HST/ACS images helps 
to examine the form of the scaling relations at faint luminos- 
ity and low surface brightness, and to investigate how the 
Coma dwarfs fit onto the well-known fundamental plane. In 
addition, we study the scatter of faint dEs about the funda- 
mental plane of giant ellipticals and examine whether it de- 
pends upon photometric and structural parameters of galax- 
ies, such as Sersic index, luminosity and colour. Moreover, 
we examine the form of the photometric plane of our sample 
dEs. 

This paper is arranged as follows, ^describes the pho- 
tometric data analysis and modelling the surface brightness 
of galaxies. Sj3]is devoted to study the scaling relations of our 
sample galaxies, such as the fundamental plane, photometric 
plane, size-surface brightness, and size-luminosity relations, 
as well as the correlation between different photometric and 
kinematic parameters. Departure of faint dEs from the FP of 
bright elliptical is discussed in ij4] In §5\ we study the vari- 
ation of the M/L ratio of the galaxies in our sample across 
their total mass. Discussion and summary are drawn in sec- 
tions ^6]and J7] 

Throughout this paper, we assume that the Coma clus- 
ter is located at a distance of 100 Mpc (z=0.0231), which 
corresponds to a distance modulus of 35.00 mag and angu- 
lar scale of 0.463 kpc arcsec -1 for h — 71, Q, m = 0.27 and 
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Table 1. The photometric and dynamical parameters of early-type galaxies in Coma cluster of which 28 galaxies have spectroscopic data from 
DEIMOS spectrograph. The velocity dispersion for 34 galaxies were extracted from MG05 and for 9 galaxies from Co09. The photometric 
parameter are derived from ACS images taken by filter F814W. is the central surface brightness of each galaxy found by the IRAF task 
ellipse and fig is the same but for galaxy model created by Galfit. R e is the effective radius of the galaxy encloses half of the total light of the 
galaxy, n and cr (column 10 and 11) are the Sersic index and velocity dispersion respectively. C (column 12) is the concentration parameter 
which is defined as C = 5 logioi^so / f2o) ■ 
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SI a = 0.73 (Carter et al. 2008). In this paper, all magnitudes 
are in the AB system. 



2 THE PHOTOMETRIC DATA 

The HST ACS Coma cluster treasury survey is a deep two- 
passband imaging survey of one of the nearest rich cluster of 
galaxies. The completed survey covers 274 arcmin 2 area of 
sky in the core and infall region of the Coma cluster. 25 fields 
were imaged by ACS Wide Field Camera with the F475W 
(g-band) and F475W (I-band) filters. Of 25 fields, 19 were 
located within 0.5 Mpc (0.3 deg) of Coma centre. For the 
purpose of this study, we used the images of Data Release 
2 (DR2), which include several improvements to the initial 
release. 



2.1 Photometry of Sample Galaxies 

Among the galaxies with DEIMOS spectroscopic data, 32 
have HST/ ACS images. Reliable velocity dispersions were 
derived for 28 of these galaxies. Four remaining galaxies 
were therefore excluded from the analysis. Moreover, ve- 
locity dispersion measurements were available for 41 more 
galaxies (34 from MG05 and 9 from Co09) in Coma for 
which ACS images at the above bands were available. The 
final sample consists of 71 galaxies with ACS images and 
determined velocity dispersion covering a luminosity range 
fromMfl w -22 to -15. 

For the purpose of this study, galaxies in both 
F814W/F475W images were extracted using SExtractor 
(Bertin & Arnouts 1996) for photometry and measuring the 
initial model independent shape parameters such as the ef- 
fective radius, R e , position angle and ellipticity. For each 
galaxy, the initial central position as well as the concentra- 
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tion parameter (i.e. C = 5 logio(rso/r2o), where rso (V20) 
is the radius within which %80 (%20) of the total galaxy 
light is collected) and the Kron radius, i.e. a characteris- 
tic radius as weighted by the light profile originally defined 
by Kron (1980), were measured using SExtractor. For each 
galaxy, we used both F814W/F475W bands to run SExtrac- 
tor in dual-image mode where F814W-band was used for ob- 
ject detection. The SExtractor input parameters are directly 
taken from Table 1 of Hammer et al. 2010. 



2.2 Galaxy Surface Brightness Profile 

Sersic function, defined as log IOOocr 1 /™, describes the struc- 
ture of most elliptical galaxies remarkably well (Kormendy, 
2009). To find the best Sersic fit to the light profile, we used 
Galfit (version 3, Peng et al. 2010). For ACS images, the 
anisotropic PSF shape depends on the location of each object 
on ACS CCD chips (WFC1 & WFC2) and was modelled by 
tinytim (Krist, 1993). To run Galfit, initial values of the R e , 
/i e , galaxy position and galaxy position angle were taken 



from SExtractor initial run, and the initial value of the Sersic 
index n was set to 3. Tests show that with the well defined 
psf of ACS data the final solution does not depend strongly 
upon the initial value of n (see also Hoyos et al. 201 1). Any 
object in the vicinity of the target galaxies was masked out. 
It is important to leave enough sky background for a reliable 
estimation of the background level, as the estimated Sersic 
index, n, is slightly sensitive to the masked regions. This 
sensitivity is higher for larger Sersic indices. Comparing the 
estimated effective radius, (R e , i.e. the radius encompass- 
ing half-light of the galaxy) and effective surface brightness, 
(( / u) e , i.e. the mean surface brightness within R e ), from Gal- 
fit and SExtractor, we found that the difference in results (the 
scatter around the line with the slope "one") is minimized 
when the size of the fitted area (galaxy and background) is 
about 2.5 times the Kron radius. In each case, the modelled 
galaxy and the Galfit residual images were inspected by eye 
to identify galaxies that are well described by Sersic model. 
The galaxies with internal spiral structure or those with poor 
fit were excluded from the analysis. 

To study the dependency of measured Sersic indices 
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on the observing wavelengths, surface brightness fitting was 
performed on both F475W/F814W images. The Sersic in- 
dices in F814W band are about 10% higher than those of 
F475W band (see Figure[TJ. To study the wavelength depen- 
dency of the FP and PHP in i ]3.2| for each band, the corre- 
sponding Sersic index and effective radius were used. 

In addition to Galfit and SExtractor, the IRAF task 
isophot in STSDAS package is used to find and compare the 
radial light profile of each galaxy and its corresponding es- 
timated Sersic model. This helps to examine the reliability 
of the fitted profiles. Moreover, the central surface bright- 
ness of galaxies, fj,o, are obtained from their estimated sur- 
face brightness profiles. Figure [2] shows the example of two 
galaxies modelled by Sersic function. GMP 3080 seems to 
have an extra component at its centre which is well modelled 
by additional Gaussian function with FWHM ~ 0.13", in 
another iteration. The importance and physical meaning of 
this excess light is discussed in £|4.2| 

The estimated photometric and kinematic parameters 
of our 71 sample galaxies are presented in TableQ] Hoyos et 
al. (2011) present a detailed comparison of effective radius 
and surface brightness with the ground-based data of Gutier- 
rez et al. (2004) and Aguerri et al. (2005). They find a good 
agreement, with a few outliers where complex structure is 
not well resolved in the ground-based data. Our derived val- 
ues are in general in good agreement with those of Hoyos et 
al. (2011) although, our derived values of the Sersic index 
are somewhat lower. This difference can be attributed to the 
fact that we allow a separate nuclear component in some of 
the fits. 



3 THE SCALING RELATIONS 

Taking all essential kinematic and photometric parameters, 
we have investigated the most well known scaling relations 
for our sample galaxies. Our sample consists of 71 dwarf 
galaxies which are fainter and less massive than the pre- 
viously studied galaxies in the Coma cluster. The distribu- 
tion of magnitude, velocity dispersion and Sersic index of 
our galaxies are represented in Figure [3] In this section, we 
also compare the scaling relations of our sample of dEs with 
other studies. 

3.1 The Fundamental Plane (FP) 

In order to fit the FP to our sample galaxies, we performed 
an orthogonal distance regression, in which the sum of the 
orthogonal residuals about the FP is minimized. Compared 
with the ordinary least-square fitting method, in which the 
residuals in (/i) e or R e are minimized, the orthogonal dis- 
tance regression is less sensitive to the outliers and is more 
robust (JFK96; La Barbera el al. 2010). During the fitting 
process, the square inverse of the error-bars of the measured 
quantities were used as weight numbers, and therefore the 
final fitting results are less affected by the values with large 
error-bars. To do the weighted orthogonal distance regres- 
sion, we used the ODRPACK software (Boggs et al. 1989, 
1992). 

The conventional projection of the FP is shown in Fig- 
ure [4] The fitting was performed on bright galaxies (i.e. 



Table 2. Coefficients of fundamental plane, FP for different lu- 
minosity cutoffs. The FP relation is represented as log(R e ) = 
A log(a) + B < n > e +C. 





14W 


A 


B 


RMS 


< - 


-20 


1.33±0.02 


0.32±0.00 


0.031 


< - 


18 


0.77±0.11 


0.17±0.03 


0.116 


< - 


-16 


0.97±0.16 


0.19±0.02 


0.150 



(/x)f 814W " < 21). In panels c and e, we have excluded the 
galaxies with poor Sersic fits. There is no significant differ- 
ence in the fitted FPs in panels (a) and (c) within la uncer- 
tainty. 

We probe the dependency of FP on the wavelength by 
performing the Sersic fit to galaxies in F475W filter. In H2A\ 
we showed that the measured Sersic parameter in red filter 
(F8 14W) is % 10 higher than the same in blue filter (F475W). 
Panel (e) in Figure [4] shows the FP for the same galaxies in 
F475W ACS (g-band) filter. The resulting FP for both red 
and blue filters (panels c & e) are compatible with each other 
within the error bars which implies no colour dependency of 
the FP of our dEs. 

For the 12 brightest galaxies (M$i4, < —20 mag) in 
our sample, the best orthogonal fit for FP is derived as 

log(R e ) = (1.33 ± 0.02) log(a) + (0.32 ± 0.00) {fi) e 
+ (-8.74 ±0.11), (2) 

with a RMS scatter of 0.031 dex in log(i? e ). R e and a 
are in kpc and km s _1 , respectively. Even though, our study 
is based on the precise measurements of the structural and 
photometric parameters from the HST/ACS imaging data, 
the coefficients of FP in this analysis for bright galaxies 
is consistent with those derived from optical (V-band) and 
near-infrared (K-band) data for 48 giant ellipticals in the 
Coma cluster (Mobasher et al. 1999: Mo99). Moreover, our 
resulting FP for bright ellipticals agrees well with other stud- 
ies of the FP. For instance, JFK96 found the coefficients of 
log(a) and {fi) e as 1.3 and 0.32, respectively. In addition, La 
Barbera et al. (2010) have investigated the FP for a sample 
of local early-type galaxies (z < 0.1) from the SDSS DR6, 
with M r < —20, in the grizYJHK wavebands. In agree- 
ment with our results, they have shown that the coefficient of 
(/i) e is independent of the waveband and equals to 0.32 by 
minimizing the orthogonal residuals about the FP. The RMS 
scatter about the FP in previous study of bright ellipticals in 
the Coma cluster are 0.074 dex in log(_R e ) (Mo99) for op- 
tical V-band, 0.08 dex (JFK96; Djorgovski & Davis 1987) 
and 0.07 dex (de Carvalho & Djorgovski 1992) for the B- 
band and 0.075 dex (Lucey et al. 1991) for the V-band, using 
other independent samples. The RMS scatter about the fitted 
FPs depends on the sample completeness and the magnitude 
range of the studied galaxies. We noticed that adding fainter 
galaxies results in larger scatters, indicating that faint dwarf 
galaxies do not lie on the FP of their brighter counterparts. 



The best coefficients of FP relation for 35 galaxies of 
our sample, brighter than M 8 i 4 = -18, are A=0.77 ± 0.11 
and B=0.17 ± 0.03 which result in a RMS scatter of 0.116 
dex in log(i? e ) corresponding to an uncertainty of 28 per 
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Figure 3. The distribution of galaxies in our sample in terms of the absolute magnitude in F814W-band, central velocity dispersion, cr, and 
Sersic parameter, n. 



cent in distances to individual galaxies. Considering all 
galaxies in our sample with Afgi4 < — 16 results in a FP 
with A=0.97 ± 0.16 and B=0.19 ± 0.02 with a RMS scatter 
of 0. 150 dex in log(R e ), which corresponds to an uncertainty 
of 35 per cent in distances to individual galaxies. The co- 
efficients of FP for faint magnitude cutoffs, Msn < — 18 
and Msi4 < —16, differs from those of equation [2] for 
-M814 < —20. This implies that, compared to bright el- 
lipticals with A/814 < —20, the faint dEs lie on an en- 
tirely different FP with larger scatter. In addition, We found 
that the fitted FPs for two sample of faint dEs with differ- 
ent faint-end cutoffs (i.e. Msi4 < — 18 and Afsi4 < —16) 
are consistent with each other within the error bars, suggest- 
ing the same FP for dwarf galaxies in the magnitude range 
-20 < A/ 8 i4 < -16 mag. In Table |2] the results of the 
best fitted FPs for different magnitude ranges are presented. 
We study the departure of faint dEs from the FP of bright 
ellipticals in S|4] 

We also studied the FP in ft-space (Bender et al. 1992). 
k parameters form an orthogonal coordinate system which 
combine the central velocity dispersion, the effective radius 
and the effective surface brightness. In this formalism, k1 
and k3 are proportional to log(M) and log(M / L) respec- 
tively, in which M is the mass of the galaxy while L is its 
luminosity. Moreover, K2 is proportional to log(M/ L){n)l 
and measures the galaxy compactness for a given mass. The 
projection of the galaxies in K-space is illustrated in Fig- 
ure[5] The left panel of Figure [5] shows that the M/L ratio 
of faint dEs does not decrease on the same linear trend as 
brighter ellipticals. The variation of the M/L ratio across 
our sample galaxies and its physical meanings are studied 
in [j5] As seen in the middle panel of Figur^U the occupied 
region of our sample dEs in the n\ — K2 diagram agrees well 
with the Figure 3 of Maraston et al. (2004). This indicates 
that dEs are generally form a linear trend perpendicular to 
the occupied region by bulges and giant elliptical galaxies. 
The different behaviour of giant and dwarf ellipticals in the 
ki — K2 space indicates that they are distinct families with 
different formation scenarios. Similar behaviour, albeit with 
smaller samples of lower-luminosity galaxies, is seen in Fig- 



ure 9 of Guzman et al. (1993) and in Figure 1 of Aguerri & 
Gonzalez-Garcfa (2009). 



3.2 The Photometric Plane (PHP) 

The Sersic index and the central velocity dispersion of galax- 
ies are correlated (see Figure [6} and hence, one is able to 
use the Sersic index instead of the velocity dispersion when- 
ever the kinematic parameters of the galaxies are not avail- 
able. Replacing the velocity dispersion, a, with the Sersic 
index, n, in FP relation (EquationQJ, the photometric plane 
(PHP) is obtained more economically than the fundamen- 
tal plane, only based on the photometric properties of the 
galaxies. Very similarly, PHP could be valuably used as a 
diagnostic tool to study the galaxy evolution and constrain 
the early-type galaxy luminosity evolution with redshift (La 
Barbera 2005: Lb05). In order to provide a local universe 
reference and to compare the relation between FP and PHP 
of the early-type galaxies, we construct the PHP for our sam- 
ple of dwarf galaxies. 

Right panels of Figure [4] show the projection of best 
PHP for our sample galaxies regardless of their brightness. 
We note that most of the outliers in top right panel have inter- 
nal spiral structure or poor Sersic fit. It appears that while the 
fundamental plane is less sensitive to detailed morphology 
of the galaxies, the photometric plane is able to effectively 
differentiate the morphological type of the galaxies. In right 
middle panel, we have eliminated the outliers of top panel in 
order to only have elliptical galaxies with good Sersic light 
profiles (hereafter "the modified sample"). This reduces the 
RMS scatter in {[i) e from 1.12 to 0.74 mag. In addition, we 
have constructed the PHP using F475W ACS (g-band) data 
in order to examine the wavelength dependency of PHP. The 
resulting PHP for both filters (panels d & f) are compatible 
with each other within the uncertainties. This is also consis- 
tent with the results of multi-band (R/I/K-band) study of the 
PHP of early-type galaxies (Lb05). Considering the follow- 
ing representation for the PHP, 
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Figure 4. Fundamental and Photometric planes (FP & PHP) in different filters, a) The dash-dotted line represents the projection of the FP fitted 
for galaxies with (yu) e < 21. Most galaxies with {fi)e > 21 tend to have higher velocity dispersion compared with that expected from the FP 
of bright galaxies. R e and a are in terms of kpc and km s^ 1 , respectively. High errors occur when measuring the velocity dispersion of fainter 
galaxies (e.g. (fx) e > 21) with very low signal-to-noise-ratio in their spectra. Therefore, in this diagram, the horizontal error bars are bigger 
than the vertical error bars, b) The dash-dotted line represents the projection of the Photometric Plane fitted for galaxies with {/i) e < 23. Most 
of the outliers are SO galaxies or galaxies with internal structures. These galaxies typically do not have very good Sersic fits. All red asterisk 
outliers belong to dSO or SB galaxies, c) Same as the panel (a) eliminating the galaxies which are not well modelled by Sersic function and have 
internal structures. The modified sample consisting of dwarf ellipticals, dEs, are represented in this panel. To avoid any confusion, the error 
bars are not shown. The typical size of error bars are the same as panel (a) d) Same as the diagram (b) for the same sample as in panel (c). The 
typical size of error bars are the same as in panel (b). e and f) Same as the panels (c) and (d) for F475W-band. All photometric properties and 
profile parameters (e.g. Sersic index) were derived from F475W images. The resulting FP and PHP in both filters, show no discrepancy within 
la error bars. The difference in coefficient "c" comes from the different magnitude of each galaxy in different pass-bands. In all panels, open 
black circles represent dwarf ellipticals from our DEIMOS/Keck observations. For red asterisks and green filled circles, the velocity dispersions 
were obtained from Cody et al. (2009) and Matkovic & Guzman (2005), respectively. 
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Figure 5. Dwarf elliptical galaxies of the Coma cluster plotted in re-space (Bender et al. 1992). All galaxies with bad Sersic fit, dSO galaxies and 
those with internal spiral structure are not plotted in this diagram. Ki, K2 and K3 are proportional to log(M), log(M / L){fj)^ and log(M/ L), 
respectively. The velocity dispersion for black open circles were measured using our DEIMOS data. For red asterisks and filled green circles, the 
velocity dispersion extracted from Co09 and MG05, respectively. In left panel, the dash-dotted line, (K3 = (0.20±0.06)ki + (— 0.09±0.17)), 
is the best fit for galaxies with {/i) e < 20, for which rei is typically greater than 2.50. The estimated linear trend for bright galaxies in our 
sample is consistent with the results of Treu et al. (2006) who found ft3 = (0.21 ± 0.02)fti + (0.19 ± 0.08) in B-band, for a large sample of 
Coma galaxies. This diagram indicates that the M/L ratios of less massive dEs are larger than those predicted by the extrapolated line to lower 
values of rei (rei <x log(M)). The dashed diagonal line in the middle panel (rei + re 2 = 8) represents the border of the "zone of avoidance" 
not populated by early-type galaxies (Burstein et al. 1997). 



Table 3. RMS scatter of the data points about the best fitted fundamental and photometric planes (see Figure [4). The scatters are calculated 
along the effective surface brightness, {fi e ) and effective radius, R e . For brighter galaxies ((/i) e < 21 in F814W-band) the scatter about 
the fundamental plane (FP) is smaller than the corresponding photometric plane (PHP). Considering all data points or only fainter galaxies 
((/^}e > 21 in F814W-band) indicates that PHP reduces the RMS scatter in both (/^ e ) and R e . 



Fitted Plane 


Filter 


Sample 


RMS scatter along {fi) e [mag] 


RMS scatter along log(R e ) [dex.] 


all data 


{n)e < 21 


(M > 21 


all data 


We < 21 


{n)e > 21 


FP(a) 


F814W 


all 


1.01 


0.35 


1.50 


0.32 


0.11 


0.47 


FP(c) 


F814W 


modified 


1.04 


0.28 


1.49 


0.31 


0.08 


0.45 


FP(e) 


F475W 


modified 


1.04 


0.28 


1.49 


0.34 


0.09 


0.49 


PHP (b) 


F814W 


all 


1.10 


1.24 


0.87 


0.38 


0.43 


0.30 


PHP (d) 


F814W 


modified 


0.72 


0.72 


0.73 


0.30 


0.30 


0.31 


PHP(f) 


F475W 


modified 


0.76 


0.60 


0.92 


0.36 


0.28 


0.43 



log{R e ) = a log{n) + /3 {u) e + 7, (3) 

Lb05 obtained a ~ 1 and j3 ~ 0.2, both in the op- 
tical and near-infrared wavebands, for a sample of galaxies 
brighter than Mi w -20 mag located in "MS 1008-1224" 
galaxy cluster at z — 0.306. Lb05 have reported the intrin- 
sic RMS dispersion in R e as 32%. Re-deriving the PHP co- 
efficients for the modified GR02 sample of galaxies in Virgo 
and Fornax clusters and using the same regression method, 
Lb05 found a and /3 as 1.14±0.15 and 0.180±0.031, re- 
spectively. This implies no redshift dependence of the PHP 
of bright elliptical galaxies. 

For galaxies brighter than Msu ~ Mr ~ — 18 mag 
in our modified sample, we found a = 1.82 ± 0.10 and 
/3 = 0.27±0.02 withaRMS scatter of 0.173 dexinlog(ii e ), 
corresponding to an uncertainty of 40 per cent in distances 
to individual galaxies. Considering all galaxies brighter than 
Mi w -16 the results are a = 2.78±0.23 and /3 = 0.39± 
0.04, with a RMS scatter of W dex in log(i? e ) or 61 per cent 
error in distances to each galaxy. 

In this study, we have few bright galaxies and there- 



fore we cannot compare our fitted PHP for bright galaxies 
with results in the literature. The measured coefficients of 
the PHP relation for our dEs differs from those of previous 
studies of the PHP of of bright elliptical galaxies (GR02; 
Lb05). We attribute this discrepancy mainly to the difference 
in mass range, most of our galaxies have a < lOOfcms -1 
whereas previous samples have a > lOOfcms -1 . The typ- 
ical Sersic indices of our galaxies are less than ~2 while 
those of GR02 and Lb05 are greater (see Figure [3} . As an- 
other possibility to explain the discrepancies, we note that 
our study is based on the HST/ACS photometric data and, 
therefore the measured photometric and structural parame- 
ters are likely to be more accurate than those measured from 
ground based observations. 



3.3 Fundamental Plane vs. Photometric Plane 

As Figure|4]shows, the deviation of galaxies from the FP and 
the PHP become larger as their luminosities, central velocity 
dispersions and Sersic indices decrease. In order to quantify 
the dispersion around the planes, the RMS deviation about 
the fitted planes are listed in Table[3] The deviations are de- 
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Figure 6. The relation between Sersic index, n, of the galaxies un- 
der study and the photometric and dynamical parameters (i.e. Mgi4, 
/irj, (A*)e, concentration parameter and cr). (jj,) e and [iq are the ef- 
fective and central surface brightness of the galaxies, respectively. 
For definition of concentration parameter refer to i]2. II All photo- 
metric parameters are extracted from F814W ACS images. For open 
circles, the velocity dispersions are derived from DEIMOS spec- 
tra in this study and for red asterisks and filled green circles they 
are extracted from Co09 and MG05 catalogues, respectively. Open 
stars belong to the galaxies with bad Sersic fit or with internal spiral 
structures. In each panel, 'a' and 'b' are respectively the slope and 
intercept of the fitted lines using orthogonal regression. CC is the 
cross-correlation coefficient of the fitting parameters. 



rived in terms of the mean effective surface brightness and 
the effective radius for three cases: (i) all data points, (ii) 
brighter galaxies with {u) e < 21 and (iii) fainter galaxies 
with {fi) e > 21. The table shows that the scatter around the 
PHP is less than the corresponding FP when we consider all 
of our modified sample galaxies (panels c and d of Figure[4j 
or when galaxies are fainter than ({fJ,) e > 21). For our mod- 
ified sample (i.e. galaxies with good Sersic fit or without 
any internal structure) the RMS scatter in {/j,) e for the FP is 
1.04 mag while it is reduced to 0.72 mag for the PHP. The 
corresponding deviation in R e is 0.31 dex for the FP and 
0.30 dex for the PHP. If we only consider the fainter galax- 
ies (i.e. (fj) e > 21), the RMS deviation in (fi) e is 1.49 mag 
for FP, while PHP decreases this value to 0.73 mag. The 
RMS scatter of brighter galaxies (i.e. (fj) e < 21) around 
the PHP increases compared to the corresponding FP. This 
implies that the brighter galaxies are still better placed on 
the FP. We conclude that faint galaxies agree better with the 
photometric plane of brighter galaxies, compared to the fun- 
damental plane. This implies that the galaxy substructures 
which are not reflected in velocity dispersion, and hence, in 
the FP relation, are better reflected in the Sersic parameter, 
and hence, in the PHP relation. In addition, the larger errors 
in velocity dispersion measurements for fainter dEs may be 
responsible for their larger scatter about the conventional FP 
compared to the PHP. 



3.4 Correlations with the Sersic index 

In Figure [6] we present the correlations between the Sersic 
index of the sample galaxies and M%ia, (fj)e, Mo, c an d 
concentration parameter. We find a linear trend between 
the Sersic index and the central surface brightness as 
log w {n) = (2.00 ± 0.35) - (0.09 ± 0.02) fi with the cor- 
relation coefficient of -W. This is entirely consistent with 
the relation between these quantities found by Graham & 
Guzman (2003). Replacing n by a we find a similar rela- 
tion (i.e. log w (a) = (3.97 ± 0.56) - (0.12 ± 0.02)^ ) 
with weaker correlation coefficient of -0.80. We ignore the 
outliers in the fitting process, which are the galaxies with 
poor Sersic fit, represented by open asterisks. Figure [6] in- 
dicates another correlation between the central velocity dis- 
persion of our sample galaxies, a, and their Sersic indices, 
n, with the correlation coefficient of 0.79. This correlation 
enables us to construct a relation a and n as logio(a) — 
(0.59 ± 0.14)k>5in(ra) - (0.77 ± 0.20) where cr is in km 
s" 1 . 

The correlation between the light concentration in 
galaxies and their Sersic index is also presented in the 
bottom panel of Figure [6] Ignoring the outliers which are 
illustrated with open asterisks, this relation is written as 
log w (n) = (0.36±0.06)C*-(0.98±0.18) with the correla- 
tion coefficient of 0.95, where C is the concentration param- 
eter and is defined as C — 5 ^oc/io(rso/)"2o)- The process of 
determining the Sersic index is model dependent while the 
concentration parameter is model independent and well de- 
termined using a simple photometric analysis. The C — n re- 
lation shows that "C" and "n" could be interchangely used. 
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Figure 7. The linear relation between (/Lt) e and log(R e ), also 
known as the Kormendy relation, for different magnitude ranges in 
F814W-band. The coefficient "p" is the slope of the linear trends: 
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Figure 8. The effective radius against the luminosity of our sample 
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best linear trend which is found using an orthogonal distance regres- 
sion (dashed line). As seen, R e and luminosity are correlated with 
the correlation coefficient of CC=-0.70. The RMS scatter about the 
fitted line is 0.17 dex. along the vertical axis, R e . 



3.5 The size-surface brightness relation 

The linear relation between the effective surface brightness, 
(/x) e , an d half-light radius, R e , of elliptical galaxies, also 
known as Kormendy relation (KR, Kormendy 1977), is rep- 
resented as 

(n) e =plog(R e )+q. (4) 

Figure [7] shows the best fit of KR for our dEs in differ- 
ent luminosity ranges, Msu < —20, — 20 < Mgi4 < —18 
and —18 < M$i4 < —16. For our brightest galaxies 
(M 8 i4 < -20), we find the KR slope as 4.27±0.18 which 
is significantly comparable to 2.43±0. 15 for a sample of gi- 
ant galaxies in Coma cluster with central velocity disper- 
sion a > 200 km s _1 (Ziegler et al. 1999). We attribute 
this discrepancy to the fact that our sample covers differ- 
ent range of size and magnitude compared with the galaxies 
studied by Ziegler et al. (1999). The typical velocity disper- 
sion of our dEs is less than 100 km s _1 and the average 
size of our dEs is ~1.5 Kpc. Moreover, the effective radii 
of our galaxies with Msi4 < —20 are less than ~4.5 Kpc, 
while all other studies of KR cover galaxies with larger sizes 
(>~ 10 Kpc) in this magnitude range (e.g. La Berbera et 
al. (2010); Ziegler et al. 1999; D'Onofrio et al. 2006). The 
imposed systematic restrictions on our sample of dEs, such 
as luminosity and size cuts, changes the geometric shape of 
the distribution of galaxies on the log(i? e )-(/^) e plane, and 
therefore results in different KR slope (Nigoche-Netro et al. 
2008). Fitting the KR to our dEs in the magnitude range 
-20 < M 8 i4 < -18 and -18 < A/ 8 i4 < -16, we found 
the KR slope as 5.23±0.37 and 5.17±0.28, respectively. In 
agreement with our results, Khosroshahi et al. (2004) have 
also found the slope of KR as ~5.2 for dwarf ellipticals of 
16 nearby galaxy groups with — 14 < Mr < —18. 

The KR is originated from the definition of the effective 
radius, R e , which relates the luminosity and effective sur- 



face brightness as L = 2tvI £ R^. Theoretically, p equals 5, 
and any difference from 5 is the results of the change in ge- 
ometric shape of the distribution of galaxies on the \og(R e )- 
(/i) e plane. Any change in magnitude range of the galaxies 
and the shape of the magnitude distribution results in differ- 
ent slope of KR (Nigoche-Netro et al. 2008). 

In agreement with Khosroshahi et al. 2004 and 
D'Onofrio et al. (2006), Figure [JJ shows that for fainter 
galaxies the log(i? e )-(/i) e linear trend is steeper than that of 
the brighter galaxies (Msi4 < —20). We also noticed that 
the slopes of KR for our galaxies in the magnitude range 
-20 < M 8 i4 < -18 and -18 < M 8 i4 < -16 are consis- 
tent within the error bars. Due to the limited R c range, we 
did not fit the KR for galaxies with — 16 < Msi4. 



3.6 The size-luminosity relation 

In Figure[8] we show the size-luminosity relation, R e oc L T , 
for our sample dEs. We found r « 0.24, which is consistent 
with the results of de Rijcke et al. (2005) who found the B- 
band radius-luminosity power-law slope between 0.28 and 
0.55 for a sample of dwarf ellipticals and dwarf spheroidal 
galaxies, r depends on the luminosity range of the galaxies 
under study. As Desroches et al. (2007) found, r systemat- 
ically varies from r ~ 0.5 at M r ~ —20 to r ~ 0.7 at 
M r « —24, which are in general in agreement with our re- 
sult if one extrapolates them to lower luminosities despite 
the structural differences between luminous and faint dEs. 



4 DEPARTURE OF GALAXIES FROM THE FP 

In order to understand the deviation of the faint dEs from 
the FP associated with the brighter galaxies, we study the 
correlations between the departure of galaxies from the FP 
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Figure 9. The deviation of the real effective surface brightness, (jj.) e , from that predicted by FP, (/i) * , against different dynamical and 
photometric parameters. Panel (A) shows the tight correlation between the deviation from the FP and the luminosity of the galaxies. Panel 
(B) represents weaker correlations in terms of n. In panel (C), the deviation is drawn in terms of the concentration parameter which is model 
independent and defined as 5 * logio(rgo/r2o)- F° r a 'l cases we used the least square method to find the best linear trend (dashed lines) by 
minimizing the residuals in App. 'a' and 'b' are slope and intercept of the fitted lines, respectively. 



with other observables (e.g magnitude, Sersic index, con- 
centration parameter, velocity dispersion, colour). In this 
study, we define the departure of galaxies from the FP as 
App — |(/i) e — (a*)J|> where (fi) e is the effective surface 
brightness of the sample galaxies and (fi)* is the expected 
effective surface brightness from the FP of bright ellipticals. 

4.1 Dependence on Galaxy Light Profile 

In Figure [9] the correlation between App and magnitude, 
Mgi4, Sersic index (n), and the concentration parameter are 
studied. The data points are binned and the error bars are 
based on the lcr RMS scatter of Afp within each bin. As 
different panels of this Figure show, A pp is anti-correlated 
with the luminosity of galaxies, their Sersic indices and their 
light concentration. Although these three parameters are cor- 
related with each other (see Figure[6]l, this graph shows that 
the maximum correlation is seen between the luminosity of 
galaxies and App with a correlation coefficient of 0.65. 

In this study, we obtain the partial correlation coef- 
ficients between App and other observables, defined as 
P cc (App, a; 0). Here, P cc measures the correlation be- 
tween Afp and one of the observables, a, in the case where 
the influence of the third variable, ft, is eliminated. We 
obtained the partial correlation coefficients between App 
and luminosity by taking out the effects of a and n as 
P CC {A FP , M 8 i 4 ; tr)=0.64 and P CC (A FP , M S14 ; n)=0.45. 
In addition, we found that Afp is very weakly correlated 
with a with the correlation coefficient of -0.30. This is con- 
sistent with the fact that a itself contributes in the FP re- 
lation and hence no strong correlation between Afp and 
a is expected. As another fact, taking out the effect of a 
(or n), when calculating the P cc (App, Mgi4', a or n), does 
not significantly change the correlation coefficient between 
App and Msu- Furthermore, the correlation coefficient be- 
tween App and n is -0.52 and after taking out the effect of 
galaxy magnitudes, we obtain P cc (App,n; Msi4,)=-0.19. 



We also obtained P cc (App, n; <r)=-0.46. These all imply 
that the departure from the FP is correlated more tightly with 
luminosity than with the Sersic index. 

The departure from the FP, App, is also correlated with 
the light concentration in the sample galaxies, C, the cross 
correlation factor of -0.56. As the concentration is expected 
to be correlated with the intrinsic luminosity of the galaxies, 
we found partial correlation to be P cc (App,C; Mgi4)=- 
0.17. This is similar to the App correlation with the Sersic 
index, as both indicate roughly the same in galaxy light pro- 
file. This can also be inferred from the tight correlation be- 
tween Sersic index and concentration parameter (see Figure 

HD>. 

We have also examined the correlation between App 
and effective radius, R e , of our sample galaxies, as a repre- 
sentative of their size. In this case, the correlation coefficient 
is obtained as -0.30. This implies that App is almost size in- 
dependent. This is also expected form the contribution of R e 
in the FP relation. 



4.2 Dependence on Galaxy Central Excess Light 

Another possible parameter which may be correlated with 
App is the central excess light relative to their best single 
Sersic model (see Figure[2]l. This central excess light (CEL), 
or nuclear star cluster, occurs in the majority of dwarf galax- 
ies in the Coma cluster (den Brok et al. 201 1), and is some- 
what bluer than the main body of the galaxy. The excess 
light is an imprint of the formation history of galaxies and its 
colour and distribution indicate that it is a result of some dis- 
sipational process. A probable explanation is that it is the re- 
sult of the last major merger, in which the gas content of the 
gas-rich progenitor falls into the core and undergoes star for- 
mation generating the excess light (e.g. HCH08, Kormendy 
et al. 2009). Numerical simulations (Makino & Hut 1997) 
show that the galaxy merger rate scales as N 2 /a 3 , where N 
is the density of galaxies, and a is velocity dispersion. Hence 
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Figure 10. Deviation from the FP against the relative central excess 
light (i.e. CEL, see the text). CEL is measured within an aperture 
of 0.4" diameter (i.e. ~0. 19 kpc). The black open diamonds show 
the galaxies in our DEIMOS sample. For red asterisks and filled 
green circles, the velocity dispersion were extracted from Co09 and 
MG05, respectively. The most deviant galaxies are labelled with 
their GMP IDs (1983). 



in Coma, with 3.5 times the density and a 1.2 times greater, 
mergers should occur 7 times more frequently than in Virgo, 
den Brok et al (201 1) conclude that the observed relationship 
between Sersic index and the strength of the colour gradient 
also points towards a history of wet mergers in the history of 
dwarf galaxies. 

Alternatively the nuclear excess could be formed from 
gas already in the dwarf galaxies as they fall into the clus- 
ter. Ram pressure stripping will remove gas from the outer 
parts and tidal interactions could drive the residual gas into 
the core, where it forms stars. In this case we would expect 
an environmental dependence of the CEL properties. Such 
dependence will be complex, as both the stripping which re- 
moves gas, and the tidal interactions driving residual gas to 
the core, will be strongest in the centre of the cluster, un- 
fortunately due to the failure of ACS in 2007, and the con- 
sequent truncation of the HST survey, the sample of dwarfs 
with high spatial resolution images outside the cluster core 
is small. 

An inspection the radial surface brightness profile of 
our sample galaxies confirmed that most part of the galax- 
ies, beyond the aperture with 0.40 arcsec diameter, are well 
modelled by Sersic function. Therefore, we measured the 
central excess light of our sample galaxies within apertures 
with diameter of 0.4 arcsec (i.e. 0.19 kpc at Coma distance), 
which were co-centred with the target galaxies. Here, we 
represent the central excess light as 

CEL = ( F ^~ F mod \ ^ (5) 

\ Fmod J aper 

where, F ga i is the flux of the target galaxy within the 
defined central aperture, and F mo d is the flux of the fitted 
Sersic model, within the same aperture. The extra flux is 



Table 4. The cross correlation coefficient of the colour of our sample 
galaxies and other photometric and dynamical parameters. Please 
refer to ^4.3l for definition of colour parameters. Figure fTTI also rep- 
resents App versus the Colour parameters. 



Correlation 


App 


Mgi4 


logio(a) 


logio(n) 


Coefficients 


(1) 


(2) 


(3) 


(4) 


g-r 


-0.48 


-0.73 


0.57 


0.72 


Am 


-0.54 


-0.78 


0.58 


0.74 


Am(04) 


-0.58 


-0.78 


0.61 


0.74 


Am(10) 


-0.59 


-0.81 


0.62 


0.78 



plotted against the deviation of galaxies from the FP (App) 
in Figure[l0] Even though the App — CEL correlation co- 
efficient is 0.12 and we did not find any clear relation be- 
tween App and CEL. In general, the most deviant galaxies 
(the labelled points) tend to have more central excess light. 
As Figure [10] shows, the relative CEL of all galaxies with 
App > 1 mag are more than 10%. 

It is likely that these galaxies have extra luminosity at 
their central regions due to star formation imposed by either 
wet mergers of gas driven to their centres by tidal interac- 
tions. The star formation activity expels out the luminous 
baryonic matter from the galaxy and changes its mass-to- 
light ratio, M/L. Any change in the M/L ratio results in 
the deviation of the galaxy from FP. We study the M/L ra- 
tio of our sample dEs in Sj5] 

4.3 Dependence on Galaxy Colour 

The location of early-type galaxies relative to the FP is a 
result of their age and stellar population properties (Graves 
et al. 2010; Graves et al. 2009, Terlevich & Forbes 2002; 
Forbes at al. 1998). Here, we study the colour- App relation, 
in order to assess the stellar population dependency of the 
faint dEs deviation from the FP of bright ellipticals. 

The colour of our dEs are calculated in 4 forms. First, 
we used the SDSS DR7 data to construct g — r for all 
of our galaxies. Second, we used our HST/ACS images 
in two F814W (I-band) and F475W (g-bang) to calculate 
Am = ?ri475 — m8i4. In addition, we estimated the galaxy 
colours within two co-centred apertures with diameters of 
0.4 and 1.0 arcsec (equal to 0.30 and 0.75 kpc), located at 
the centre of each galaxy (hereafter, Am (04) and Am (10), 
respectively). In table[4] we present the correlation factor for 
each pair of defined colours of our sample galaxies and other 
photometric and dynamical parameters (i.e. Mai4, cr, n and 
App). 

Table [4] shows that App is tighter correlated with the 
colour of galaxies compared to their central excess light (see 
i]4.2| >. In all cases, Am(10) displays stronger correlation 
compared to all other defined colour parameters. Moreover, 
Am (10) is stronger correlated with the galaxy magnitudes, 
Mgi4, implying that fainter galaxies are bluer. A strong cor- 
relation also exists between Am(10) and logio(n), which 
indicates that bluer galaxies have smaller Sersic indices and 
therefore are less concentrated. 

In Figure [TT] the deviation of galaxies from the FP is 
plotted against different colour indicators. As seen, the bluer 
galaxies have larger App values. In agreement with our re- 
sults, Bernardi et al. (2003) have also found a correlation be- 
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Figure 11. Deviation of galaxies from the FP (F814W-band), App, against their colours. The colour of galaxies in panel (A) were calculated 
based on SDSS DR7. In panel (B), colours were derived from our analysis on ACS F814W/F475W images. Panel C & D are the same as panel 
(B) except for colours which were derived within the apertures of 0.4" and 1" diameter at the centre of galaxies. At the bottom of each panel, 
all data points were binned. Each bin contains the same number of galaxies. For all cases we used the least square method to find the best linear 
trend (dashed lines) by minimizing the residuals in Afp. 'a' is the slope and 'b' is the intercept of the fitted lines. CC is the cross-correlation 
coefficient of the fitting parameters. This diagram shows that the bluer galaxies are more deviant from the FP specifically when considering the 
colour of galaxies at their 0.4" or 1.0" central regions. The most deviant galaxies are labelled with their GMP IDs. 



tween the the residuals to the FP and colour of their sample 
of bright early- type galaxies. Figure 3 of their paper indi- 
cates that bluer galaxies display greater App. 

This implies that the deviation of dEs from the FP is 
correlated to their stellar population. The supernovae activi- 
ties, as a consequence of star formation in bluer galaxies, has 
swept away the luminous baryonic matter from these galax- 
ies. Thus the supernova driven winds truncate star formation 
and modify the M/L ratio of the bluer galaxies, increas- 
ing their scatter about the FP. This scenario also explains 
the location of the faintest galaxies of our sample on the L- 
a diagram (see Paper I), indicating that these galaxies have 
larger velocity dispersions compared to the trend of brighter 
dEs. We noted that, these galaxies are bluer than the other 
galaxies in our sample and display larger deviation from the 
FP. The star formation have ejected the luminous matter of 
these galaxies and shifted them to the faint end of the L-cr 
diagram, while their total mass (mostly the dark matter con- 
tent), and hence, their internal velocity dispersions remain 
the same. 



5 THE M/L RATIOS 

The behaviour of M/L ratio may account for FP tilt (e.g. 
ZGZ06; Cappellari et al. 2006; Grave & Faber 2010) and/or 
deviation of galaxies from the FP (e.g. Reda, Forbes & Hau 
2005). Based on the virial theorem, the dynamical mass, M, 
within the effective radius, R e can be derived as, 

M S c4^, (6) 



where a is the galaxy central velocity dispersion and G 
is the gravitational constant. Assuming that our sample is ho- 
mologous, c is almost constant for our galaxies and depends 
on the profiles of the dark and luminous matters. We adopt 
c = 5 (Cappelari et al. 2006); although our conclusion is 
independent of the exact value of c as long as it is uniformly 
applied to all galaxies. In this study, we assumed that the ab- 
solute magnitudes of the Sun (in AB system) in F475W and 
F814W bands are 5.14 and 4.57 mag, respectively. 

In Figure [12] we show the variation of the dynami- 
cal mass-to-light ratio, M/L, as function of galaxy dynam- 
ical mass, M. For galaxies brighter than Mgi4 = —17.5 
mag, the correlation between M/L and M (the dotted line) 
is obtained as M/L oc M 0,23±0 04 , in agreement with 
the results of other studies for bright ellipticals. As an ex- 
ample, Mo99 found M/L oc M° 18±0 ' 01 in K-band and 
M/L oc M°' 23±0 ' 01 in optical V-band, for galaxies in the 
mass range of 10 10 -10 13 M Q . In addition, studying a sam- 
ple of 9000 local early-type galaxies of SDSS in the red- 
shift range 0.01 < z <0.3, Bernardi et al. (2003) found 
M/L oc M~°' 2 . 

In our sample of dEs, the best linear relation between 
logio(u) and {fi) e is found to be logio(a) — (4.94 ± 
0.75) - (0.16 ± 0.02) (/x) e with the correlation coefficient 
of -0.66. This relation is translated to I e oc o - 2 - 5±0 - 3 ! where 
I e is effective surface brightness in flux units, in F814W- 
band. Using the results of Binggeli, Sandage & Tarenghi 
(1984), we expect J e oc a -2 ' 5 for elliptical galaxies. Using 
the virial theorem for spherical systems, a 2 oc GM e /R e , 
and the relation between luminosity and effective surface 
brightness, L — 2nI e R 2 , the M/L ratio is obtained as 
M/L oc a 2 /y/LTL (Co09). Substituting our derived Faber- 
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Figure 12. The mass-to-light ratio (M/L) of our sample galaxies 
in terms of their dynamical mass. Panels (A) and (B) are based on 
ACS F814W-band and F475W-band photometric data, respectively. 
As seen, the slopes of (M / L)-M lines do not depend on the two 
passbands (F814W and F475W). Open circles represent the galaxies 
with DEIMOS data for which their a were measured in paper I. Red 
asterisk and green filled circles belong to dEs from Co09 and MG05, 
respectively. The dotted line is the best linear fit for brighter galaxies 
with A/814 > —17.5. 'a' and 'b' are the slope and the intercept of 
the dotted lines, respectively. In panel (A), the typical error bars of 
M/L are illustrated. The mass, M, and M/L ratios are in solar 
units. 



Jackson relation (Paper I), L oc 1.99 ± 0.14, and the de- 
pendency of I e upon cr, we obtain M/L oc (j- - 25 * - 33 x 

L -0.13±Q.17 K M -0.15±0.22 for galaxies in our sam . 

pie. Our result is more consistent with Co09 who found the 



M dyn /L oc M" 



for a sample with almost the same 



luminosity range as ours. 

As seen in panels (A) & (B) of Figure [T2l for Msi4 > 
— 17.5, some galaxies have larger M/L ratios with respect 
to the faintward extrapolation of the linear trend of brighter 
galaxies, while few faint galaxies follow the M/L — M 
relation of bright galaxies. The most deviant galaxies in 
M/L — M diagram (see Figure \\2\ are the bluer galax- 
ies, too. This suggests that the formation of the most de- 
viant dwarf galaxies and those which are following the trend 
of bright ellipticals can be explained by different scenar- 
ios. Nevertheless, we still need more data point to study 
the scatter of faint dEs about the trends of brighter ellipti- 
cals and to examine different formation mechanisms at the 
faint regime. In agreement with our results (see Figure IT2t. 
studying the dwarf galaxies with — 16 < Mv < —12 in 
the core of Perseus cluster, Penny et al. (2009) found that 
fainter dwarfs have larger M/L ratios. In addition, Geha et 



al. (2002) showed that fainter galaxies have larger M/L ra- 
tios compared to their brighter counterparts 

We have also investigated the waveband dependency 
of (M/L)-M relation, to explore whether it is governed 
by the change in metalicity of galaxies and/or their stellar 
population. Comparing panels (A) and (B), which are based 
on ACS F814W and F475W-band photometry, we find the 
same behaviour of M/L ratio in different wavebands. We 
noted that, the (M / ' L)-M relation of galaxies fainter than 
M$i4 = —17.5 is also independent of the passband. There- 
fore, the stellar population is not the only parameter respon- 
sible for the change of M/L ratios (Mo99). 

As a conclusion, we found that the M/L ratio is not 
constant over all of our sample dEs and varies with the mass 
and the luminosity of the galaxies. The variation in M/L ra- 
tio is responsible for the deviation of our fainter dEs from the 
FP. Since the fainter galaxies in our sample are bluer than the 
other galaxies and have larger M/L ratios, we attribute their 
deviation from the FP to their recent star formation activities 
(see also $43}. 



6 DISCUSSION 

In this paper, we studied the fundamental and photometric 
planes of a sample of 71 dEs in the core of Coma cluster, 
the nearest massive elliptical-rich cluster down to luminos- 
ity of A/814 < -15.3. Taking advantage of the DEIMOS 
high resolution spectrograph, which enables us to measure 
the internal velocity dispersion of dwarf ellipticals, and high 
resolution imaging of HST/ACS which allows an accurate 
surface brightness modelling, we were able to extend the FP 
of galaxies to ~ 1 magnitude fainter than the previous stud- 
ies. 

We obtained the FP for a subsample of 12 
galaxies brighter than Msi4 = — 20 as R e oc 

a 1.33±0. 02^-0. 80±0. 01 which {f , consistent with the pre . 

vious studies of bright galaxies in Coma (JFK96; Mo99). 
Studying the FP of 141 early-type galaxies in the Shapley 
super cluster at z=0.049, Gargiulo et al. (2009: Ga09) found 
that the FP follows the relation R e oc ct 135 (7)-° 81 for 
galaxies with a > 100 km s _1 and Mr < —18.7. When 
including all galaxies in their sample, including low-mass 
galaxies down to a ~50 km s~ , Ga09 found a shallower 
exponent for a. The FP relation of our dEs displays even 
shallower exponents for a and {L) e than in Ga09, due to fur- 
ther extension to fainter galaxies. 

In Figures [4] and QT| the faintest galaxies and the most 
deviant data points are based on our DEIMOS observations. 
One could argue that our a measurements might be overes- 
timated. In paper I, we have considered all source of un- 
certainties in measuring the velocity dispersions, namely 
the statistical errors, template mismatch and other system- 
atic uncertainties. We have also checked the sensitivity of 
our measured velocity dispersions to different kind of stel- 
lar templates. Moreover, using different set of mixed stel- 
lar templates, covering a vast range of spectral types, when 
measuring the velocity dispersions, confirmed our measure- 
ments. The presented uncertainty in the measurement in- 
cludes all such sources and therefore we can rule out the 
overestimation argument. 
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In this study, the spectroscopy was carried out using 
the slits with the size of 0.7 arcsec covering the body and 
the central region of the galaxies. Due to the size of slits, 
for galaxies with net rotation, the measured velocity disper- 
sion differs from the central velocity dispersion. The kinetic 
energy per unit mass (KE) for a spheroid is (e.g. Busarello, 
Lango &Feoli 1992): 

KE = \i v ) = \ V ?ot + ^oo ~ !<r, 2 n , (7) 

where V ro t is the rotation velocity, (To is the central ve- 
locity dispersion and a m is the measured velocity disper- 
sion. Thus we have 

(£=L)a~ 1 + \(Yl2Ly. (8) 

(To 6 (To 

The observed values of V ro t/<To range from 0.01 to 
~1.0 for galaxies with — 14 < Mr < — 26 (Davies et al. 
1983; Pedraz et al. 2002; Simien & Prugniel 2002; Geha 
et al. 2003; van Zee et al. 2004), implying that the differ- 
ence between a m and oo is less than ~15% of the central 
velocity dispersion. This kind of uncertainty is smaller than 
the total estimated uncertainty in measuring the velocity dis- 
persions of faint dEs (a < 50 km s -1 ). Based on the re- 
lation for an oblate, isotropic galaxy flattened by rotation 
(Davies et al. 1983), the mean isophotal ellipticity of four 
galaxies in our sample, GMP2780, GMP2563, GMP3119 
and GMP3141 are less than 0.1, corresponding to V ro t/cro 
~ 0.3. Using equation[8] the uncertainty in the measured ve- 
locity dispersion due to the rotation would be ~2%, which 
is significantly smaller than the uncertainties due to template 
mismatch and other systematics discussed in paper I. The 
uncertainty due to the rotation is ~6% for GMP2655 and 
GMP2808, for which the mean ellipticity is ~0.25. Thus, ro- 
tation can not be solely responsible for the deviation of faint 
dEs from the faintward extrapolation of the Faber- Jackson 
relation and the FP of bright ellipticals. 

The non-universality of the FP is claimed by ZGZ06, 
where they suggest that the spheroids ranging from dEs to 
the galaxy clusters lie on a curved surface, in (<r,i? e ,(/x) e ) 
space. Describing the M / L ratio as a parabolic function of 
the internal velocity dispersion, ZGZ06 could explain the 
change in the coefficients of the FP for different spherical 
systems. They predicted that for a ~ 50 km s~ , M/L is 
almost a flat function of a, with a large scatter, which tends 
to increase for smaller values of a. In agreement with Co09 
and ZGZ06, we found M/L oc ,j-o- 2 5±o-33 for aR galax _ 
ies in our sample. In addition, we found larger M/L ratios 
for faintest dEs in our sample with respect to the (M / L)-M 
trend of brighter galaxies (see Figure [72t . 

A mixture of formation mechanisms has been consid- 
ered to model dE formation. In the "wind model", dEs are 
primordial objects that lost their gas in a supernova-driven 
galactic wind (DS86; Yoshii & Arimoto 1987: YA87; Chiosi 
& Carraro 2002: CC02). In contrast to massive ellipticals, 
dwarf ellipticals have very long star formation histories. 
Since giant galaxies hold on very strongly their gas content, 
their gas is almost completely converted into stars in a sin- 
gle burst. On the other hand, supernova activities disperse 
the gas and switch off further star formation in dwarf galax- 



ies. As the gas cools, it sinks back in and again undergoes 
a new phase of star formation followed by supernova explo- 
sions which disperse the gas, and this cycle can be repeated 
many times (Bender & Nieto 1990; de Rijcke et al. 2005). 

An alternative formation path for dwarf galaxies in- 
volves stripping of larger galaxies by gravitational and gas 
dynamical processes. As galaxies fall in to the cluster they 
are stripped of much of their gas by ram pressure effects 
(Gunn & Gott 1972; van Zee, Skillman & Haynes 2004), 
this process can be seen at work in both Ha and Ultra- Violet 
imaging of the Coma cluster (Yagi et al. 2010; Smith et al. 
2010). Stars and dark matter can then be removed by in- 
teractions with more massive galaxies and the cluster po- 
tential, termed "harassment" (Moore et al. 1996; Moore et 
al 1998; Mayer et al. 2001; Mastropietro et al. 2005) or by 
tidal encounters between galaxies of near-equal mass (Rich- 
stone 1976; Aguerri & Gonzalez-Garcfa 2009). Both harass- 
ment and tidal encounters can create bars (Mastropietro et al. 
2005; Aguerri & Gonzalez-Garcfa 2009) which could play a 
part in driving gas not stripped by the ram pressure of the in- 
tergalactic medium inwards, triggering star formation. The 
resulting galaxy is likely to have a higher central surface 
brightness than its progenitor and would probably resemble 
a nucleated dE which rotates quite rapidly and still displays 
some memory of its former state. Although harassment tends 
to increase internal velocity dispersion, it is unable to disrupt 
rotational motion and therefore is not obviously reconciled 
with the low rotational velocities found by Simien & Prug- 
niel (2002), Geha et al. (2002, 2003) and MG05. 

Another possible way to explain the formation of dwarf 
galaxies is presented by Kroupa (1998) and Due et al. 
(2004). The merging process of two gas-rich disk galaxies 
results in many massive star clusters which coalesce to pro- 
duce a small number of dwarf galaxies. These tidally formed 
dwarf galaxies contain negligible amounts of dark matter 
and hence cannot explain a significant fraction of the dE 
population. 

Our findings appear to be most consistent with the 
"wind model". In the wind-stripping model, the shallow po- 
tential well of dwarf galaxies, with the velocity dispersion of 
less than 100 km s _1 does not allow intense star formation 
without stripping most of the gas (Schaeffer & Silk 1988). 
We found the effective radius-luminosity and Faber-Jackson 
relation as R e oc L 2i and L oc cr 2 in F814W-band, re- 
spectively. These two relations are in agreement with the cal- 
culations of de Rijcke et al. (2005) based on the winds mod- 
els of YA87 and CC02. They showed that, although, YA87 
models do not predict a tight L — R e relation, CC02 models 
can reproduce it very nicely. 

Using the accurate measurement of the M / L ratios for 
the SAURON galaxy sample (Bacon et al. 2001), the FP tilt 
of early-type galaxies can be attributed to the variation of 
M / L ratio which is governed by the stellar population and 
dark matter properties (Cappellari et al. 2006). Moreover, in 
agreement with Cappellari et al. (2006), considering the re- 
sulting M I L ratios based on different star formation models 
and mass profiles, Allanson et al. (2009) and Grillo & Gobat 
(2010) concluded that non-homology is not a major driver 
of the FP tilt while the mass sensitivity of the M/L ratio 
is the primary source of FP behaviour in different mass and 
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luminosity regimes. We find that the fainter dEs are bluer, 
and scatter more about the FP, than the brighter ellipticals. 
Star formation activity, together with the wind model, can 
explain how the bluer dwarfs at a given mass have different 
M/L ratios and thus follow different scaling relations. 



7 SUMMARY 

The results of this study can be summarized as follows: 

(i) The FP of bright ellipticals with M 8 i 4 < -20 in 
our sample is fully consistent with the previous studies of 
the Coma bright galaxies. We found that the scatter about 
the FP depends at the faint-end luminosity cutoff is signifi- 
cantly increased. We noticed that the FP tilt of faint dEs with 
Msi4 >-20 differs from that of giant ellipticals. 

(ii) The FP is not colour dependent which implies that 
the stellar population plays a insignificant role in the scat- 
ter about the FP and changing its tilt. We studied the de- 
parture of faint dEs from the FP of bright ellipticals, Afp, 
which is found to be correlated with the luminosity and 
the galaxy brightness profile. Less concentrated faint dEs, 
smaller Sersic indices, show larger deviation from the FP. 
We also found a relation between Afp and galaxies central 
excess light as an imprint of their formation history. In addi- 
tion, we found a correlation between Afp and the colour of 
galaxies which indicates that the bluer galaxies in our sam- 
ple dEs display larger deviation from the FP of the bright 
ellipticals. On the other hand, we have already shown that 
in the L-er relation, discussed in paper I, fainter dEs seem to 
be more massive (e.g. dynamical mass estimated using the 
internal velocity dispersion) than predicted from their lumi- 
nosity. Although tidal effects can be important in removing 
luminous matter from the galaxy, we believe that supernova 
driven winds are more probably the dominant process. 

(iii) Replacing the central velocity dispersion with the 
Sersic index in the FP relation, we studied the photometric 
plane, PHP, of our sample galaxies. The differences in the 
coefficients of PHP relation for our dEs and those of previ- 
ous studies in local elliptical galaxies (GR02), mainly arise 
from the fainter luminosity coverage of our sample. We no- 
ticed that most of the outliers in the projected PHP diagram 
have internal morphological features or poor Sersic fit. Our 
conclusion is that, while the fundamental plane is less sen- 
sitive to detailed morphology of galaxies, the photometric 
plane is able to effectively differentiate between different 
morphologies. In other words, the scatter about the FP in- 
creases as we extend the sample to fainter dwarf galaxies. 
In the PHP, however, the scatter is driven by asymmetry and 
unsmooth galaxy brightness distribution. 

(iv) We studied the correlations between the Sersic in- 
dex, n, and other photometric or kinematic parameters (i.e. 
Mgi4, (^i} e , pio, <? and the concentration). The light concen- 
tration and Sersic index are strongly correlated with the cor- 
relation coefficient of 0.95. The process of determining the 
Sersic index, n, is model dependent. On the other hand, the 
concentration parameter is model independent and is deter- 
mined using a simple photometric analysis. Therefore, based 
on the n-concentration correlation, one is able to economi- 
cally replace the Sersic index, n with the concentration pa- 
rameter when processing the scaling relations. 



The best radius-luminosity relation for our sample dEs 
is obtained as R e oc L 0,24 which is consistent with the 
results of de Rijcke et al. (2005) who reported a B-band 
radius-luminosity power-law slope between 0.28 and 0.55 
for a sample of dwarf ellipticals and dwarf spheroidal galax- 
ies. 

(v) We found the Kormendy Relation (KR) slopes as 
5.23±0.37 and 5.17±0.28 for galaxies in the luminosity 
ranges of -20 < M 8 i 4 < -18 and -18 < M 814 < -16, 
respectively. This is in agreement with the results of Khos- 
roshahi et al. (2004) who found the KR slope as ~5.2 for 
dEs of 16 nearby galaxy groups with — 18 < Ms < —14. 
Furthermore, in agreement with Khosroshahi et al. (2004), 
we found the steeper linear trends for faint dEs compared to 
that of the brighter galaxies. 

(vi) We obtained the mass-to-light ratio as M/L oc 

a -0.25±0.33 x L -0.13±0.17 ^ M -0.15±0.22 for our sample 

dEs indicating that the M/L ratio is almost a flat function of 
mass, luminosity and velocity dispersion. Our finding agrees 
well with that of Co09, who reported a M/L oc M 09±0 06 
for a sample with nearly the same luminosity range. Galax- 
ies with Mgi4 > —17.5 seems to have higher M/L ratios 
with respect to the linear extrapolation of the same associ- 
ated with brighter galaxies. The variation in M/L ratio is 
responsible for the deviation of the fainter dEs from the FP. 
Since, the fainter galaxies in our sample are bluer than the 
other galaxies and have higher M/L ratios, we attribute their 
deviation from the FP to their recent star formation activi- 
ties. 
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